Fermiology, Orbital order, Orbital Fluctuation and Cooper Pairing in Iron-based 

Superconductors 
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We address two important issues that arise in recent studies of iron-based superconductivity. (1) 
Why are the T c of Aa;Fe2-i/Se2 and the single unit cell FeSe on SrTiC>3 so high despite both only 
have electron pockets? (2) What (if any) are the effects of orbital order and orbital fluctuation on 
the Cooper pairing. Our conclusions are summarized in the third paragraph of the paper. 



The discovery of A x Fe 2 -ySe 2 [l} (T* 



48K, s-wave pairing in A a .Fc 2 - !; Sc2 and FeSe/STO. 



under pressure [2]) and single unit cell FeSe on 
SrTi0 3 (FeSe/STO)[3] (T™ aa; =65K, determined by 
angle-resolved photoemission spectroscopy (ARPES)[4]), 
stirred up a new wave of excitement in iron-based super- 
conductors (FeSCs) research. In ARPES studies it is 
found, at ambient pressure, both systems have no hole 
pocket [4, 5]. Because it is often perceived that the scat- 
tering between the electron and hole pockets are impor- 
tant for both antifcrromagnctism and Cooper pairing[6], 
this becomes an issue. 

On a different front, recently many experimental ev- 
idences point to the fact that FeSCs have a tendency 
to become electronically "nematic" [7-9, 11, 12]. For ex- 
ample, through magnetic torque measurement Ref.[ll] 
reported a phase diagram for BaFe2(Asi_ a: P a ;)2 where 
the superconducting dome is enclosed by a non-magnetic 
electronic nematic phase. In addition, Rcf.[12] reported a 
divergent "nematic susceptibility" in Ba(Fei_ a; Co 2; )2As2 
close to the x value at which T™ a:E occurs. In addition to 
these, an ARPES experiment by Yi et al.[9] established 
the tie between electronic nematicity and the d xz , d yz or- 
bital ordering. These experiments naturally raise the 
question: what role (if any) do orbital order or orbital 
fluctuation play in Cooper pairing? 

The purpose of this paper is to address the above two 
questions. Our conclusions are summarized as follows. 
(1) Hole pockets introduce frustration in Cooper pairing 
(a concept we shall discuss later). To a large extent 
this is due to the existence of band vorticity around 
the hole fermi surface (see later). Removing the hole 
pockets releases pairing frustration, experimentally it 
is found this docs not weaken the antifcrromagnctic 
(AFM) correlation[13]. This makes Cooper pairing in 
A a; Fe2-; / Se2 and FeSe/STO stronger, hence higher T c . 
Of course in FeSe/STO substrate screening can further 
enhance T c [14]. (2) Orbital fluctuation has negligible 
effect on Cooper pairing while static orbital order can 
have large effect. (3) AFM fluctuation can still be the 
primary cause of Cooper pairing, but AFM is most likely 
due to local correlation not fermi surface nesting. (4) 
The interpocket hybridization tends to favor in phase 



The above conclusions arc reached by using an effective 
Hamiltonian approach. In this approach we write down 
a low energy Hamiltonian to capture the system's ten- 
dencies toward (1) stripe antiferromagnetism[15-17] (2) 
superconducting pairing[18, 19], and (3) d xz /d yz orbital 
ordering[7-12]. The above tendencies are not only estab- 
lished experimentally, but also found theoretically[6, 20- 
23]. To a great extent, the derivation of the effective 
Hamiltonian has been achieved by the functional renor- 
malization group (FRG) study in Ref.[25]. It's conclu- 
sion has been checked by variational Monte-Carlo (VMC) 
calculation[26]. The last statement is significant because 
FeSCs are by no means weakly coupled systems [27]. 

Our effective Hamiltonian is given by: 

H c ff = ^ EkQ^kacr + ^ Fkaa n kaanqi3T + 

k,a,fr k,a,<rq,/3,T 

(i) 

ij i 

The first term is the bandstructure (a, f3 = band in- 
dices, er, r = spin indices). The second term is a 
multi-band version of the Fermi liquid interaction. In 



the third term S. t = lJ2a c i,a,s°ss'Ci, a , s >, where a = 
d z 2,d xz ,d yz ,d X y,d x 2_ y 2 is the orbital index. The last 
term of Eq. (1), with V > 0, describes the tendency to- 
ward the d xz /d yz orbital ordering[9, 23, 24]. We wrote 
the last two terms of Eq. (1) in real space, but they should 
be understood as been projected to the band eigen bases 
that lie in a thin shell around the fermi surface. 

Eq. (1) describes several competing (or "intertwined") 
instabilities, a hallmark of strongly correlated systems. 
Of course solving the effective Hamiltonian is a difficult 
problem. In the following we shall assume Cooper pairing 
is the winning instability, and our goal is simply to deter- 
mine which pairing symmetry is favored the most. Some 
technical points: (1) In Eq. (1) we allowed the magnetic 
interaction to extend over arbitrary neighbors. However 
in Ref. [25] it has been shown if one retains only the first 
(Ji) and second neighbor (J2) interaction, the effective 
Hamiltonian already qualitatively captures the numerical 
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FIG. 1: (a) A caricature of the fermi surface of Fe-pnictides. 
Here blue and red mark the hole and electron pockets, re- 
spectively, (b) The two regimes of pairing. The value of the 
critical angle (marked by the red line) is approximately 0.37T. 
In the region marked "F" pairing is frustrated, (c) The typ- 
ical gap function in region "I" of panel (b). There is hatch 
size is proportional to the magnitude and the color indicate 
the sign (blue: plus and red: minus). Only the first quadrant 
of the unfolded Brillouin zone is shown. 



functional renormalization group results for the k depen- 
dence of the AFM and SC order parameters. Since the 
purpose of this paper is to elucidate qualitative physics 
rather than providing quantitative predictions, we shall 
truncate to only J\ and J2P8]. (2) Once the fermi 
surface is fixed, the second term of Eq. (1) has no effect 
on Cooper pairing. Hence we shall drop it in the sub- 
sequent discussions. Finally on the semantics: when the 
various interaction in Eq. (1) fail to drive long range or- 
der, we call them "fluctuations". For example the J\, J2 
terms will be termed "magnetic fluctuation" . 
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FIG. 2: The gap function in region "F" of Fig. 1(b). In the 
presence of the hole pockets (left panels in (a) and (b)) the gap 
function has s-wave symmetry and the sign of the gap function 
on the hole pockets changes as 9 increases from (Fig. 2(a)) 
to 0.3 7r (Fig. 2(b)). After removing the hole pockets (right 
panels in (a) and (b)) the symmetry becomes d-wave. Only 
the first quadrant of the unfolded Brillouin zone is shown. 
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FIG. 3: The band vorticity associated with the hole pockets, 
(a) The band structure near F. (b) The pseudo spin winding 
(see text).(c) The region with pairing frustration substantially 
shrinks after the removal of the band vorticity of the hole 
Fermis surfaces. The values of the two critical angles (marked 
by the red lines) are approximately 0.15-7T and 0.287T. Only the 
first quadrant of the unfolded Brillouin zone is shown. 



Results for Fe-pnictide: A plot of the typical fermi 
surface for the Fe-pnictide systems is shown in Fig. 1(a), 
where both electron pocket (marked red) and hole-pocket 
(marked blue) are present. First we study the effect of 
magnetic interaction on pairing. We control the ratio 
between the AF J\ and J2 by introducing an angle 9 

J 1 =Jcos9, J 2 = Jsw9, 0<9<n/2. (2) 

We then project the effective Hamiltonian onto the sin- 
glet pairing channel to construct the following pairing 
"matrix" 

M{a,<j>;M') = U{a,4>;W) ^fjy (3) 

VFp.rW) 

Here U(a,4>; {},(/>') is the effective singlet pairing inter- 
action on the fermi surface (a labels fermi pockets and 
4> is the angle around them), kpp and VFp,r are the 
magnitudes of the fermi wavevector and the projection 
of the fermi velocity along the radial direction, respec- 
tively. The leading gap function is the eigenfunction of 
M(a, </>; P, 4>') with the most negative eigenvalue. In case 
M (a, </>;/?, </>') has degenerate eigenvalues non-quadratic 
interaction terms between different superconducting or- 
der parameters are necessary to determine the pairing 
symmetry [29, 30]. In this paper we shall just focus on 
the quadratic instability. Depending on the value of 9, 
there exists two pairing regimes marked by "I" and "F" 



in Fig. 1(b). The representative gap functions in region 
I is shown in Fig. 1(c), it is s±. Cooper pairing in this 
region is non-frustrated in the sense that if we remove 
the hole pockets and recalculate the gap function on the 
electron pockets, the same sign structure is obtained. 

In contrast, in region F pairing is frustrated. In 
Fig. 2(a,b) we show the gap function with and without 
the hole pockets in this region. Generally speaking the 
pairing symmetry is s-wavc with hole pockets and d-wave 
without. The fact that after removing the hole pock- 
ets the symmetry of the gap function changes is remi- 
niscent to the change of spin alignment from coplanar 
to collinear after removing a spin to unfrustrate the AF 
Heisenberg interaction on a triangle. In the literature [31] 
it is noted that there is a near degeneracy between the 
s± and d x 2_ y 2 pairing symmetry for values of the interac- 
tion parameters yielding significant pairing strength. We 
claim this degeneracy is due to the pairing frustration 
discussed above. 

It turns out the pairing frustration is, to a large ex- 
tent, caused by the band topology of the hole pockets. 
In Fig. 3(a) we show the dispersion of the hole bands 
near the center of the Brillouin zone. The double degen- 
eracy at r is protected by the point group symmetry. A 
consequence of such degeneracy is the existence of "band 
vorticity" at T [32, 33]. To see that we first note the 
orbital content of the hole band eigenfunctions are pre- 
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FIG. 4: (a) Typical gap function in region "I" of Fig. 3(c), 
there is no pairing frustration, (b) Typical gap function in 
region "II" of Fig. 3(c), there is no pairing frustration, (c) 
Typical gap function in region "F" of Fig. 3(c). There is 
pairing frustration. The left and right panels are the gap 
functions with and without the hole pockets. Only the first 
quadrant of the unfolded Brillouin zone is shown. 



dominantly d xz and d yz . Let's use a pseudo spin 1/2 to 
represent these two orbital states, e.g., d xz — > t z = 1/2 
and d yz — > t z = —1/2. Using the band cigcnfunctions we 
compute the expectation value of the pseudospin. The re- 
sult lies in the x-z plane and the direction winds around 
the fermi surface (with vorticity 2) as shown in Fig. 3(b). 

After we remove the hole band vorticity (by modify- 
ing the band wavefunction while maintaining the band 
dispersion) the frustrated pairing region substantially 
shrinks, as shown in Fig. 3(c). The typical gap function in 
region "I", "II" and "F" are shown in Fig. 4(a,b,c). From 
this we conclude, the band vorticity associated with the 
hole pockets is an important cause of the pairing frus- 
tration. It is important to note, however, removing the 
band vorticity docs not completely eliminate pairing frus- 
tration. What's needed is the removal of the hole pockets 
entirely! 

This motivates us to make the following conjecture: 
one reason that the absence of hole pockets raises T c in 
A x Fe2- y Se2 and FeSe/STO is because the pairing frus- 
tration has been removed. However high T c also requires 
strong AFM fluctuation (at least in the magnetic pairing 
scenario). If the AFM correlation is due to the nest- 
ing between the electron and hole fermi surfaces [6], one 
would predict the weakening of AFM fluctuation (hence 
the weakening of pairing). However the recent neutron 
scattcring[13] experiment performed on "semiconduct- 
ing" Ka;Fe2-ySe2 compounds neighboring the supercon- 
ducting ones has revealed strong stripy AFM long range 
order peaking at the same wave vector as that observed 
in the parent pnictide compounds. This result suggests 
that it is better to view the AFM correlation in FeSC as 
due to local correlation physics, similar to the cuprates, 
instead of as coming from fermi surface nesting. More 
importantly it tells us the AFM fluctuation, needed for 
Cooper pairing, is strong in Kj;Fe2_j,Se2 despite the ab- 
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FIG. 5: (a) The gap function that is suppressed by the orbital 
fluctuation. (b,c) In the presence of static orbital ordering the 
frustration region also shrinks compared to Fig. 1(b). The 
phase diagram is qualitatively similar to Fig. 3(c) except the 
"F" region is slightly wider. The gap function in region "I" 
and "II" are shown in (b) and (c).Only the first quadrant of 
the unfolded Brillouin zone is shown. 



sence of the hole pockets. 

Motivated by Ref.[12] we next ask what is the effect of 
orbital fluctuation on Cooper pairing. To mimic orbital 
fluctuation we turn on the V term in Eq. (1). To our sur- 
prise, for all values of V this has no effect on pairing. To 
understand this result let us switch off Ji and J2 and look 
at the eigenfunctions of Eq. (3) when there is only the 
orbital interaction. It turns out in this case except one 
positive eigenvalue all other eigenvalues of M are zero. In 
Fig. 5(a) we illustrate the eigenfunction associated with 
the positive eigenvalue. Thus the orbital interaction lifts 
up the pairing eigenvalue associated the d xy symmetry 
but leaves all other symmetry channels unaffected. In 
particular, this lifting does not affect the leading pairing 
symmetries, namely s and d x i_ y i in Fig. 1 and Fig. 2. 
The fact that V introduces only a single positive pairing 
eigenvalue is due to the on-site nature of the orbital in- 
teraction. Had wc allowed for further neighbor orbital in- 
teraction, e.g., Vq Y<i ni, xz n lyz + Vi n i,xznj, yz there 
will be more than one positive eigenvalues. For example 
if we choose V1/V0 = 0.5 there are five positive eigenval- 
ues (all other eigenvalues are still zero). Again, none of 
these five eigenfunctions has the symmetry of the lead- 
ing pairing channel in Fig. 1 and Fig. 2. For all the cases 
we studied this is always true. Thus we conclude orbital 
fluctuation has no effect on pairing. 

Because Ref. [11] claimed the superconducting dome 
in BaFe2(Asi_ a: P 2 ;)2 is embedded in a non-magnetic 
nematic phase, we next study the effect of static 
orbital order on pairing. In Fig. 5(b,c) we show the 
anisotropic gap function in the presence of a SE = 40 
meV energy splitting between the xz and yz orbitals, 



namely, SEj^ii 7 



H,yz 



). If this anisotropy is 



observed by ARPES it will constitute a strong evidence 
for coexistence of superconductivity and nematic order. 

Results for the Fe-chalcognide: The fermiology of 
superconducting A :c Fe2- !/ Se2 and FeSe/STO are rather 
similar - there are only electron pockets. As mentioned 
earlier due to the neutron scattering result of Ref. [13] 
there is every reason to believe the same magnetic in- 
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FIG. 6: Magnetic fluctuation induced pairing symmetry of 
systems with only electron Fermi surfaces, (a) Two different 
pairing regimes. The value of the critical angle (marked by 
the red line) is approximately 0.27-7T. The gap function in 
region I (b) and II (c). Only the first quadrant of the unfolded 
Brillouin zone is shown. 
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FIG. 7: (a) The hybridization-split electron pockets, (b) The 
gap function in the presence of hybridization. In panel (b) 
the splitting between the two electron fermi surfaces has been 
enhanced to increase clarity. The value of the critical angle 
(marked by the red line) is approximately 0.1057T. 



teraction exists in the effective Hamiltonian (at least for 
K x Fe2-; / Se2 ). Using bandstructures with only electron 
pockcts[14, 34] our result for the gap function under dif- 
ferent J\ and Ji ratio is shown in Fig. 6. Depending on 
whether J\ or Ji dominates, the pairing symmetry is ei- 
ther s wave or fully gapped d-wave, consistent with the 
findings in Ref.[35]. Interestingly despite the fact that 
AFM effective interaction is characteristic of a repulsive 
system, the in-phase s-wave pairing is found for a large 
region of parameter space. 

Because nodeless d-wave occupies a significant region 
of parameter space in Fig. 6(a)), it is important to in- 
vestigate the effect of the electron pocket hybridization 
on the gap function. This hybridization exists (even for 
k z = 0) when the the z — > — z glide plane symmetry 
is broken. For A a: Fe2-ySe2 this occurs near the sample 
surface (which is what ARPES probes). For FeSc/STO 
there is no glide plane symmetry at all. It was first 
pointed out in Ref. [36] that when such hybridization is 
sufficiently strong it can modify the nodeless d-wave pair- 
ing state into nodal d-wave pairing. 

In Fig. 7(a) we choose the hybridization matrix ele- 
ments so that electron pockets is only slightly split and 
remain more or less circular. This choice is motivated 
by the fact that ARPES found rounded electron pocket 
without discernable splitting for both A x Fe2- y Se2[5] and 
FeSe/STO[4] . 

In the presence of the electron pocket splitting, the 
phase diagram and gap function are shown in Fig. 7(b). 
The hybridization has clearly increased the range of sta- 
bility for the s-wave pairing. Indeed there is a region 
of the phase diagram where the nodeless d-wave pairing 
has turned into a nodal one. However, the dominant pair- 
ing symmetry in Fig. 7(b) is the in phase s-wave rather 
than the out-of-phase s-wave proposed in Ref. [36, 37]. 
This result is consistent with the recent ARPES result 
which rules out d-wave pairing based on the mapping of 
the gap function on the electron pockets center around 
k = (0,0,71-) [38]. 

Finally we believe the substrate, namely SrTiC>3, 



plays an important role in raising the T c of FeSe/STO 
and further stabilize s-wave pairing. This is through 
the screening of the long range part of the Coulomb 
interaction[14]. As far as orbital fluctuation and orbital 
ordering are concerned, the results are very similar to 
those for systems with both electron and hole pockets. 

Conclusion: We have used an effective Hamiltonian 
to address two important modern issues raised by the 
high T c A x Fe2- y Se2 and FeSe/STO superconductors. 
This effective Hamiltonian is constructed under the guid- 
ance of (1) FRG calculation [25], (2) VMC calculation[26] 
and (3) experiments [7-9, 11, 12, 15-19]. The issues are 
(1) the effect of fermiology on high temperature super- 
conductivity and (2) the effects of orbital ordering and 
orbital fluctuations on Cooper pairing. Our conclusions 
are summarized in paragraph three of the introduction. 
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